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In an arbitrary axisymmetric stationary spacetime, we determine the expression for the tangen- 
tial velocity of test objects following a circular stable geodesic motion in the equatorial plane, as 
function of the metric coefficients. Next, we impose the condition, observed in large samples of 
disks galaxies, that the magnitude of such tangential velocity be radii independent, obtaining a con- 
straint equation among the metric coefficients, and thus arriving to an iff condition: The tangential 
velocity of test particles is radii independent iff the metric coefficients satisfied the mentioned con- 
straint equation. Furthermore, for the static case, the constraint equation can be easily integrated, 
leaving the spacetime at the equatorial plane essentially with only one independent metric coeffi- 
cient. With the geometry thus fixed, we compute the Einstein tensor and equate it to and arbitrary 
stress energy tensor, in order to determined the type of energy-matter which could produce such 
a geometry. Within an approximation, we deduce a constraint equation among the components of 
the stress energy tensor. We test in that constraint equation several well known types of matter, 
which have been proposed as dark matter candidates and are able to point for possible right ones. 
Finally, we also present the spherically symmetric static case and apply the mentioned procedure 
to perfect fluid stress energy tensor, recovering the Newtonian result as well as the one obtained in 
the axisymmetric case. We also present arguments on the need to use GR to study types of matter 
different than the dust one. 

PACS numbers: 95.35.+d, 95.35.G 



I. INTRODUCTION 



One of the most important achievements of the present cosmology is doubtless the complete identifications and 
the accounting of the amounts of the different types of matter and energy which are present in the Universe (for an 
excellent review see for example Q). Essentially, the present components of the Universe are composed by matter 
and vacuum energy Oo = Qm + [§■ Furthermore, there is a very good evidence that the Universe is flat. That 
evidence comes either from the theory, where the most accepted model for the early Universe is inflation, as well as 
from observational data, which implies that S7o = 1 ± 0.12 (see for example ||). The mass of the galaxy clusters is 
perhaps the most reliable way for determining the matter component £Im- Observations indicate that f?M ~ 0.3 [Q, 
however, the main visible components of f^M, baryons, neutrinos, form a very small fraction of £Im- Observations 
indicate that stars and dust (baryons) represents something like 5% of the whole matter of the Universe. In other 
words, Qm ~ ^fc + Ojjm ~ 0-05 + £Idm> where Qdm represents the dark matter part of the matter contributions 
which has to have a value of £Idm ~ 0.25. Notice that the given value of the amount of baryonic matter is in 
concordance with the limits imposed by nucleosynthesis (see for example A greater amount of baryonic matter 
could change the predicted values of primordial H and 4 He in the standard model of cosmology which coincide very 
well with observations. 

The existence of dark matter in the Universe has been firmly established by astronomical observations at very 
different length-scales, ranging from the local galaxies to clusters of galaxies. The standard way to notice this need for 
dark matter comes within the framework of mechanics: A large fraction of the mass needed to produce the observed 
dynamical effects in all these very different systems, is not seen. This puzzle has stimulated the exploration of several 
proposals, and very imaginative explanations have been put forward, from exotic matter like supersymmetric particles 
P| to non relativistic modifications of Newtonian dynamics |(| and non-linear general relativistic theories M, Above 
this, it is believed that this dark matter is such that interacts very weakly with ordinary matter, which makes it 
very difficult to detect by other means other than by their gravitational effects on the baryonic matter as it is well 
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stablished for the cold dark matter scenario |8|,|9|] . The bottom line is that one of the most important components of 
the density of the Universe, the dark matter, still eludes us, and the question remains open: Which is the nature of 
the dark matter component? 

At the galactic scale, the problem is clearly posed: The measurements of rotation curves (RC) in galaxies show 
that the coplanar orbital motion of gas in the outer parts of these galaxies keeps a more or less constant velocity 
up to several luminous radii (lQJLlj] . The discrepancy arises when one applies the usual Newtonian dynamics to the 
observed luminous matter and gas, since then, the circular velocity should decrease as we move outwards. The most 
widely accepted explanation is that there exists a spherical halo of dark matter, its nature being unknown, which 
surrounds the galaxy and account for the missing mass needed to produce the flat behavior of the RC. 

The main goal of this work is to study the dark matter problem in spiral galaxies, using a fully relativistic approach, 
continuing with the work started by Matos and Guzman Jl2| , where they made a preliminary dynamical analysis in the 
context of spiral galaxies. We apply a deductive method: Starting from a reasonable general space-time, we deduce, 
in terms of arbitrary metric coefficients the expression for the tangential velocity of test objects following a circular 
stable geodesic motion in the equatorial plane. We impose next the condition, observed in hundreds of galaxies, 
that such tangential velocity be radii independent, and obtain a constraint equation among the metric coefficients, 
Eq. (p5|) . Arriving in this way to an iff condition: In an axisymmetric static space-time, the tangential velocity of test 
particles at the equatorial plane is radii independent iff the metric coefficients satisfied Eq. (p5|) . Furthermore, for the 
static case, the constraint equation can be solved, leaving the space-time essentially with only one independent metric 
coefficient, thus determining very narrowly the type of space-time which can have a geometry such that the tangential 
velocity of the geodesies of test particles moving in stable circular equatorial orbits, be radii independent. With the 
geometry thus fixed, we compute the Einstein tensor and equate it to and arbitrary stress energy tensor, in order 
to determined the type of energy-matter which could produce such a geometry. We are able to deduce a constraint 
equation among the components of the stress energy tensor, Eq. (pf7|) . Using that constraint equation, Eq. (pf7|) , we 
test several well known types of matter, which have been proposed as dark matter candidates and are able to impose 
further restrictions on most of them as possible candidates for the dark matter in the region where it is observed the 
mentioned behavior on the tangential velocity. It is stimulating that the gravitational physics has been developed 
to such degree in which we can actually follow Sherlock Holmes' maxim: "...when you have excluded the impossible, 
whatever remains, however improbable, must be the truth" |]l3| . 

The work is composed as follows: In section 2 we determine the geometry of axisymmetric static space times allowing 
a tangential velocity with a magnitude independent of the radius. In section 3, we work with the Einstein equations 
for the geometry thus determined and for a general stress energy tensor, obtaining a constraint for the components of 
such stress energy tensor, and we test four types of matter into that constraint, being able to restrict most of them, 
among which is the perfect fluid. Finally, in the conclusion we discuss our results and propose a coupling of types 
of matter, which could be the one we are looking for in order to determine the nature of the dark matter. We also 
present this analysis for the spherical static space time in appendix ^ and for the axisymmetric stationary case in 
appendix fv| 



II. FORM OF THE LINE ELEMENT 

In this section we study the conditions which the flatness of the tangential velocity of the RC, imposes on the 
metric coefficients. We want to stress the fact that the results presented in this section are independent of type of 
energy-matter tensor present in the space-time and curving it. It is a purely geometric analysis. 

As mentioned above, observational data show that the galaxies are composed by almost 90% of dark matter. Thus 
we can suppose that luminous matter does not contribute in a determining way to the total energy density of the 
galaxy, at least in the region where the flatness of the RC is observed. Consequently, we consider that the dark matter 
will be the main contributor to the dynamics, and we will treat the observed luminous matter as a test fluid, that's it, 
in this approximation we will neglect the contribution of the luminous matter to the curvature, i. e., to the dynamics. 
Also, it is reasonable to suppose that the halo of dark matter is symmetric with respect to the rotation axis of the 
galaxy, thus we take the symmetry of the space-time as axially symmetric. Furthermore, the observations allow us 
to take the space-time as stationary as well. Thus, the most general reasonable space-time which we can study is an 
axisymmetric stationary one. The line-element of such space-time, given in the Papapetrou form is: 
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ds 2 = -e 2 *{dt + Ludip) 2 + e" 2 ^ [e 2l {dp 2 + dz 2 ) + fx 2 d(f 2 ], (1) 
where ip, and p, are functions of (p, z). 

We will derive the geodesic equations in the equatorial plane, that is for z = z = 0, where dot stands for the 
derivative with respect to the proper time, r. Then, we will study the constrains imposed for circular geodesies on the 
energy and angular momentum of particles in such orbits, and obtain an expression for the tangential velocity for the 
particles moving along those geodesies, described in terms of the metric coefficients. Finally, we impose the condition 
that such tangential velocity be radii independent and derive the restriction that then has to be satisfied among the 
gravitational coefficients. In this section we derive the expressions for a static space time and in the Appendix M we 
do so for the stationary one. We split the presentation in this way not only for the sake of clarity, but from the fact 
that the static case is also quite realistic, considering that the observed velocity of the stars orbiting a galaxy in the 
region of interest is quite non relativistic (of the order of 230 Km per s), thus we can infer that the space-time is not 
very rapidly rotating. 

The Lagrangian for a test particle traveling on the static space time (cu = 0) described by (|l|) is given by: 

2C = -e 2 n 2 + e~ 2 V 7 (P 2 + i 2 ) + ^ (2) 
thus, the associated canonical momenta, p x a = J^-, are: 

Pt = -E = -e 2 ^i, (3) 

Plp = L = tfe- 2 *<p, (4) 

P P = e- 2 ^p, (5) 

p z = e- 2 ^h, (6) 

where E, and L, are constants of motion for each geodesic, a fact which comes from the symmetries of the space-time 
analyzed. As there is no explicit dependence on time, t, the Hamiltonian, TL — p a x a — C, is another conserved 
quantity, which we normalized to be equal to minus one half for time-like geodesies. Also, we restrict the motion to 
be at the equatorial plane, thus z — 0. In this way, we obtain the following equation for the radial geodesic motion: 

^2 _ e 2(V-7) L<p-1]=0. (7) 

In order to have stable circular motion, which is the motion we are interested in, we have to satisfy three conditions: 

i) p = 0, and 

ii) ^M _ o, where V(p) = -e 2 ^'^ [Et-Lip-1]. 

iii) - ' \extr > 0, in order to have a minimum. 

With these conditions, from Eq. ([?]), we obtain a set of two equations constraining the motion to be circular extrema 
in the equatorial plane: 

Ei - Lip - 1 = 0, 
[Et-L<p-1]) =0. (8) 

^,From Eq.(^|), we can express t, and (p in terms of E,L, and the metric coefficients as 

i = e~ 2 ^E, (9) 

P = —L. (10) 
M 

Using these equations in the constraints ones and recalling that E and L are constants for each circular orbit, after 
some rearranging, we arrive to the following equations: 

p 2 e- 2 ^'(l-e^E 2 ) + L 2 = Q, (11) 

-( e - 2 %£ 2 +(^J) i 2 = 0, (12) 
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where the subindex stands for derivative with respect to p. Solving for E and L, we obtain: 



E = e 



\ ^~^P 



L = lie"* 



2V P 



(13) 



The second derivative of the potential V(p) evaluated at the extreme, in this case means evaluate at the values of 
E and L which constraint the motion to be circular and extrema, is given by: 



pp\extr — fj. 



pp 



iPp 



(14) 



We can now obtain an expression for the angular velocity of a test particle, J7, moving in a circular motion in the 
orbital plane, in terms of the metric coefficients, recalling that 



dt t 



thus, using Eqs.(lO), and (p_3[), in this last equation for the angular velocity, we obtain that: 

P 2V 



Q = 



(15) 



(16) 



Finally, in order to express the tangential velocity of the test particles in circular motion in the equatorial plane, 
in terms of the metric coefficients, following Chandrasekhar jl4|], we rewrite the line element given in Eq.(|l]) as: 



ds 



2 _ _2 V 



dt 2 + e" 2 * p 2 dtp 2 + e- 2 ^(dp 2 + dz 2 ), 



(17) 



thus, in terms of the proper time, dr 2 = —ds 2 , we have that 



dr 2 = e 2 ^dt 2 [1- e" 4 V (^j + 



2 / j \ 2 N 

dz 
~db 



(18) 



from which we can write that 



1 2%h Q z r 1 2i 

l = e v u II — V I, 



(19) 



where u° — 4^ is the usual time component of the four velocity, and a definition of the spatial velocity, v 2 , comes out 
naturally in this way. 



2 -4ii 2 U1 P 

v = e p, i — 



( dp 



\ dt 



( dp 



' dz 

+ 1 Tt 



(20) 



This spatial velocity is the 3-velocity of a particle measured with respect to an orthonormal reference system (see 
section 52 of [H), thus has components: 



V 2 =V M 2 +V (P) 2 +V (Z) 2 . 



^From these last two expressions we obtain for the tp— component the spatial velocity: 



(21) 
(22) 
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and substituting £1 from Eq.(16), we finally obtain an expression for the tangential velocity of a test particle in stable 
circular motion, in terms of the metric coefficients of the general line element given by Eq.(|l|), such tangential velocity 
has the form: 

(23) 

It was our goal to obtain this expression for the tangential velocity for a general axisymmetric static space time, 
and to be able to describe it in terms of the metric coefficients alone, because now we can impose conditions on this 
tangential velocity, and deduce a constraint equation among the metric coefficients, which has to be satisfied in order 
to fulfill the condition imposed on the velocity. In particular, the tangential velocity for circular trajectories in each 
orbit is constant, that is p — 0, thus = with v^ 1 a constant, representing the value of the velocity, from 
Eq. (|23|), we have that: 

+ (24 ) 

Finally, with respect to the z-motion, considering that at the equatorial plane not only z = 0, but also that z = 0, 
that is that the forces above and below the plane cancel out, from the geodesic z-equation, using Eq.(p4|), we obtain 
that this relation among the metric coefficients must hold for the derivatives with respect to z as well: 

7 = (25) 

In this way, we arrive to an if and only if condition: If Eqs. (|24| , p5|) is satisfied, then the tangential velocity of 
circular stable equatorial orbits is constant. Furthermore, if the t ang ential velocity of circular stable equatorial orbits 
is constant, then the metric coefficients have to satisfy Eqs. (|24| , |25| ). Notice then that if the function ip and fj, arc 
related by 

= (26) 
Mo 

with I = const, we obtain that this a necessary and sufficient condition for the velocity vj^ to be the same for two 
orbits at different radii at the equatorial plane, provided that I = (vj^) 2 / (l + (v c ^) 2 ). 

We call your attention to the remarkable fact that the metric coefficient 7 does not play any role in this analysis, the 
motion analyzed is determined only by the other two metric coefficients, which now are related by this last equation, 
Eq. (p5|) , thus leaving the problems in terms of only one metric coefficient. Actually this absence of 7 will be clear 
in the next section, where with the field equations we will see that it is determined in terms of the other metric 
coefficients and some components of the matter presented in the space-time, implying that it is not and independent 
function. 



Thus, in order to have tangential velocities of equatorial objects circling the galaxy, and whose magnitude is radii 
independent, the form of the line element in the equatorial plane has to be 

d f = -(jL\ > d t? + (jL) 2 [e 2 i dp 2 + (i 2 d<p 2 }. (27) 
VMo/ \[i J 

Notice that this type of space time definitely can not be asymptotically flat. Neither it has the form of a space time 
related with a central black hole. What can be said is that this line element describes the region where the tangential 
velocity of the test particles is constant all over that region, and that it has to be joined in the interior and in the 
exterior regions with other types of space times if one wishes to have a central black hole, and that the influence of 
the middle region ends at some distance and thus has an asymptotically flat external region. 

Taking into account the constraint between ip and \x given by Eq. ( p6| ) , the energy, angular momentum, the rotational 
velocity , and the second derivative of the potential have the final expressions: 
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E = S&L, (28) 
L= )p , (29) 



n = ^—(-^) , (so) 



being Z + = 1 + (u c (¥,) ) 2 and i_ = 1 - (v c (v) ) 2 . 

Notice that the second derivative of the potential at the extreme is always positive, thus the circular equatorial 
curves with constant tangential velocity are stable. 

Before going to the field equations, we think it useful to present our derivations applied to the Schwarzschild case, 
thus testing the expressions while recovering the well known results. Starting with the line elements in spherical 
coordinates: 

ds 2 = - (l - — ) dt 2 + — \ 1 dr 2 + r 2 (de 2 + sin 2 6dtf), (32) 



we perform the coordinate transformation 



v /^T^ 2 "+-^fL= + Af, fl = taa- 1 (-) J (33) 



> 2 4- z 2 z 



(the inverse transformation is p = R sin#, z = R cos 9, with R — | (r — M + \/r 2 — 2 Mr)), to obtain the line element 
in the Papapetrou form, Eq.(|), with e 2 ^ = ) , e 2 ~< = (l - T r0^ ) ) 2 ,/x = p (l - ttt&t) = 0. 



The horizon in this coordinates is located at \J p 2 + z 2 = 4p . 

Restricting the expressions to the equatorial plane, z = 0, we have that ip\ z =o = In ((p — 4£) / (p + 4f)) , A*| z= o 
p (1 — Af 2 /4p 2 ), thus from Eqs.(13), we obtain, for the energy and angular momentum: 



(p-f) 2 T ( . / M 



(p + f) v /p 2 -2Mp+M^' " V 2; y p(p2 _2Mp+^)' 
for the angular velocity test particles, from Eq.([i6|): 



= \/Mp"« I 1 + ^ ) , (35) 



2p 

which gives us the Kepler law.. For the tangential velocity, from Eq.(|2^) we obtain 

with the known dependence as the inverse of the square root of the distance. Transforming back to spherical coordi- 
nates, it can be seen that our expressions agree with the usual ones, see for instance For the second derivative 
of the potential, it is useful to write p in terms of the horizon radius, as p = n 4£, with n a number, and from Eq.(|l4|) 
we obtain: 

V pp \ eX tr- M 2 n(n + 1) 3 ( n 2_ 4rl + 1) - V<) 



G 



This second derivative is positive down to p = (5 + 2V6)4^-, which marks the last stable orbit, and corresponds to 
the known result of r — 6M in spherical coordinates. 

In this way, we are confident on our expressions and can proceed to study the field equations. 

A last remark about the geometric analysis. Recall that the observations are based on measurements of the red 
shift, not on the tangential velocity directly. If the space time is flat, the two quantities are proportional. But we are 
now working in curved space times, so we have to see if that proportionality is still valid. Following |16[] , we use the 
fact that the frequency of a photon is given by v — u a p a , with u a the four velocity of the object and p a the photon 
momentum, we have that the red shift, z, is given by 

z=l-^, (38) 

L*rec 

thus, for an object orbiting the galactic center in the equatorial plane at a distance p from the center, with tangential 
velocity v^' and emitting a photon with frequency vq, and for an observer located at rest at infinity, that is far away 
from the emission, detecting the photon with a frequency Voo, it can be shown that the red shift is given by: 



* = i- { ' + viv)) ,meL. (39) 

Vl - v^) 2 V 9tt(oo) 

Now, for the observed velocities, we have that v <C 1, i. e. they are much less than the speed of light, and we have 
to suppose that far away from the observed galaxy its gravitational influence ends, otherwise we could not detect the 
tangential velocity, we would then be moving along with the observed object! Thus we can take gtt(oo) = — 1 and we 
have that at first order in the velocity, with g u (p) = — = — 1 — ip — we obtain: 

z = -{v {ip) +ip + •••)■ (40) 

But we have computed that in the case analyzed ip — / (l + w^) 2 ) \n(p/p ) ~ 2>) 2 (/i-l)/(/i+l), thus 

ip ~ yte) , and we conclude that z ~ v^ v \ In this way, we see that the radii independence of the value of the measured 
red shift can be related with the radii independence of the value of the tangential velocity, which is the fact that has 
been studied in this work. 



III. FIELD EQUATIONS 

Now we are in a position where we can test any type of matter-energy to determine whether or not it produces 
a curvature in the space time such that the motion of the test particles can be circular stable and be such that the 
tangential velocity of those particles is constant for a large radial region in the equatorial plane. 

We obtain the general form of the Einstein tensor for the axisymmetric static space time described by Eq. ([|), 
with u> = , and equate it to an arbitrary stress energy tensor. After some manipulations we conclude that the field 
equations are a set of two equations involving the metric coefficients ip , and p: 

p{ip PP + 1>,z) +P P 4>p + Pz 4> Z =^P [e- 2 ^-^ {e'^T tt + ^ T w ) + T pp + T„], (41) 

P 

Ppp + Pzz =Snp [T pp + T zz ] . (42) 
There are also two first order equations for the other metric coefficient 7: 

IpPp ~lzPz~ P {^p 2 - Ipz 2 ) + Pzz = 8 7T pT P p, (43) 

IpPz +"fz P P - 2 pippipz - Ppz = 8 7T pT pz , (44) 

and finally, the field equations give us another equation for the second derivatives of 7 which thus is redundant, this 
equation is: 

1pp + 1zz + {^p) 2 + {^z) 2 =8tt — T w . (45) 

p 
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The analysis presented in the last section is exact and the relations between the metric coefficients and their first 
derivatives must be satisfied at the equatorial plane in order to describe the observed motion. Using the Einstein's 
equations, we need the second derivatives of those metric coefficients. Thus, we have to make the approximation that 
the relations obtained among them, holds as well in a region close to the equatorial plane. Within this approximation, 
from Eqs.(|24|, |25| , p6| ), it can be obtained the following expression: 

MVVp + + Vv + ^ = (vc iv) ) 2 / (l + («c M ) 2 ) (/ip P + Hzz)- (46) 



Thus, with this last relation, from the Einstein's equations Eqs. (|4l|,|42|) we obtain a constraint among the stress energy 
tensor components which, within the approximation made on the validity of Eqs.(|24|, ^5 26) out of the equatorial plane, 
has to be satisfied by any type of matter in order to have constant tangential velocities: 

[c ' \T pp + T zz ) = e- 2 ^(e-^T tt + ^T w ), (47) 



1 + (««>>)V v pp zz > V M 

We finish our analysis by testing several types of matter, described by their respective stress energy tensor, to see 
whether or not they are able to deform the geometry of the space time in such a way that the tangential velocity of 
the equatorial rotational objects be constant, that is, that they satisfy Eq.(E^). 



A. Vacuum Fields 



We start with the vacuum solutions, with T M1/ = 0. In this case Eq.(47) is trivially satisfied, thus we proceed 
to analyze the Einstein equations directly. From Eq.([l2|), the easiest solution implies p — p. Eq.(E0]), is a Laplace 
equation for ip and 'imposing the condition that at the equatorial plane the flat curve condition, Eq.(E5h, be satisfied, 
as well as the one of symmetry with respect to the galactic plane, we obtain that ip = I lnp and the other Einstein 
equation then imply 7 = I 2 lnp, with I — (vj-^) 2 / (l + (v c ^) 2 ). In this way, we obtained an exact vacuum solution 
for the Einstein equations, which produces that the test particles circling at the equatorial plane behave in agreement 
with the observations: 

ds 2 = -p 21 dt 2 + p- 2l (p 2 ' 2 (dp 2 + dz 2 ) + p 2 d V 2 ). (48) 

The central object is string-like. Observations show that cosmic strings object are very unlikely to exist, neverthe- 
less, this is an example of objects which could produce the observed motion of test particles, and in which the density 
does not go as r~ 2 , because it is vacuum. Thus, such a behavior on the density is a sufficient but not a necessary 
condition for the flatness of the rotational curves. 



B. Perfect Fluid 



For the perfect fluid, T pv = (d + p) u p u v + g^ v p, with d the density of the fluid and p its pressure. In this case we 
are thinking on a "dark fluid", which could be composed of planetoids or WIMPS or MACHOS, which are not seen 
but it is thought that they could be there affecting the geometry in the way needed in order to have the observed 
behavior in the tangential velocities of the luminous matter. Taking this dark fluid as static, the four velocity is given 
by u a = (u°, 0,0,0) with, for the line element given by Eq.(0) with u = 0, u° = E , and L = 0. Thus u = -E, 
and from u a u a = —1, we conclude that E = . The stress energy tensor for the dark fluid then has the form: 

T tt = dE 2 = de 2 ^ (49) 
T pp - T zz = e- 2 ^-T)p, (50) 
T vv = p 2 e" 2 V (51) 

Substituting in Eq.(^), we obtain that in the equatorial plane, in order to satisfy the observed behavior on the 
tangential velocities, the "dark fluid" has to satisfy: 

/l - (v !»1 2 \ 

-Hrrte^>= (d+P) < (52) 
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Thus, we obtain an equation of state for the "dark fluid" particles at the equatorial plane: 



p = - ^{vM? d> (53) 

which, compared to the equation of state for a perfect fluid, p = cod, implies that — 1 < lu < — |, for vj^' between 
the speed of light and zero. This result is quite remarkable. It coincides with the type of equation of state derived 
within the Quintessence model |]l7|-p0f at the cosmological level, and now we obtain similar results at the galactic 
level. This sort of matter has been called exotic matter ]2l| and studied in several contexts |2^] . Our result points 
to the fact that the Dark Matter actually could be exotic. We want to stress that due to the approximation taken 
for the behavior of the metric coefficients off the galactic plane, we are not excluding the possibility that the dark 
fluid be composed of baryonic usual matter, actually from the Newtonian approach, we know that regular matter can 
produce the observed motion. This fact is not reproduced in the present analysis, due to our approximation, what 
we certainly can conclude is that exotic type of matter also can produce such 'observed motion. 

In order to recover the Newtonian case, where we know that the dust type fluid does work as the dark matter, 
we have to analyze the spherical case, which is introduced in Appendix M. For the stress energy tensor we again 
take the static perfect fluid, thus the four velocity of the test particle reads — (it , 0,0,0), with u° — i. For the 
spherically symmetric metric, Eq. (|68|) , we know that i is associated to a conserved quantity, the energy, and is given 
by: i = -g^y, and from the normalization of the four velocity, u a u a = — 1, we get that E = \J B(r). Thus, considering 
those spacetimes for which the tangential velocity (v c ) of test particles in circular orbits is radii independent, we use 
Eq.(f76j), obtaining that u° = r - ^' , thus, u a = —T/Bor( Vo ^ . In this way, we get that the non zero components 
of the static spherically symmetric perfect fluid are: 

T tt = dB r 2(v ^, 
T rr = pA(r), 

T ee =pr 2 . (54) 



Substituting these expressions in the Einstein equations Eq.(|77|), it turns out that the equations can be completely 
solved, yielding: 

A{T) " ifT^Tij- (55) 

6Trar z Ayr) 

,= 1+ 8 2( ;f.:f (r) . (57) 

8 7r r 1 Ayr) 

with b = 2 (1 + 2 (v c ) 2 — v c 4 ), a = 1 + (v c ) 2 , and Ao an integration constant. 
Taking the particular case for the integration constant Aq = 0, we get 

Mr) = \ (58) 



(vcfi} - 

27 



d= V V , 2 2 \ (59) 
2tt br z 



Airbr 2 



(60) 



In this way, we obtain the particular solution where the particles move in the observed way, in a space-time 
with a deficit angle, g rr = constant, and in which the density goes as d ~ (v c ) 2 /r 2 , and the pressure goes as 
p ~ (i> c ) 4 /r 2 , thus it is a dust like solution. Furthermore, for an equation of state p — ujd, this dust like solution 
implies lu — (v c ) 2 / ^2 (1 — (v c ) 2 /2)j , which is between and 1, for the tangential velocity between and 1, and thus 

is a perfectly well known fluid type. In this way we recover the dust hypotheses within our approach, and we clearly 
see that the' fact that we did not recover this case within the axi-symmctric analysis, was due to our approximation 
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outside the equatorial plane. 

For the general case when the integration constant Aq is non zero, taking an equation of state as before p = u>d, 
we obtain that the to is a function of r and again, as in the axial case, it is negative, so we are dealing with the exotic 
type of perfect fluid obtained in the axi-symmetric case. 

This is a good moment to discuss the question on the need to use GR, even though the gravitational field is weak. 
In the Newtonian description it is well known that the space-time can be described as 

ds 2 = -(1 + 2<S>)dt 2 + (1 - 2$)dr 2 + r 2 dfl 2 

with $ = <I>(r) the Newtonian gravitational potential. For the spherically symmetric case we obtain that 

-gtt = r 2 ^) 2 = e 2 ^ 2 ln « = 1 + 2 (v c f ln(r) + • ■ ■ 

from here we determine that <!> = (v c ) 2 ln(r). In this Newtonian approximation the complete set of 10 Einstein 
equations reduces to one equation, the usual Poisson equation V 2 $ = AnGd. Collecting this last result, it is obtain 
that the density, d, goes as d ~ (v c ) /r 2 , which is the expression for the dark matter density known from the 
astronomers' work. 

However, notice that in this approximation it can not be said anything else about the matter producing the observed 
motion. The Newtonian approximation fixes the matter to be dust-perfect-fluid-like type. This is the usual way of 
reasoning: it is suppose a priori that the dark matter is a completely Newtonian dust and at the end of the day 
one arrives to a consistent description of the dark matter determining only the shape of the Newtonian gravitational 
potential. 

In this work we are proceeding in a different way. We are using Einstein equations backwards; we do not make any 
assumptions in the type of matter nor do approximations. From the observations on the motion of the test particles, 
we determine the geometry and then, by means of the Einstein's equations we obtain constrains of the type of matter. 
We have shown that we do recover the Newtonian result, but also it is clear that this is a very particular case for a 
very specific type of matter. In the general reasoning, we do not fix neither the type of matter nor the equation of 
state, we let the equations themselves to do that obtaining more general results. 

To end this argumentation we recall the reader that there is three conditions that have to be fulfill in order to 
reach the Newtonian limit; 1) the speeds of the study particles most be most less than the speed of light, 2) the 
gravitational field must be weak, and 3) the pressures associated with the matter study must to be must smaller than 
the corresponding density. It is this last condition which is taking a priori in the usual analysis and it is not satisfy 
in the general case as we have shown, thus justify the need of GR in order to be able to consider any type of matter. 

As a last remark about these results on the type of perfect fluid is that the Big-Bang nucleosynthesis imposes very 
strong constraints to the percentage of the baryonic matter to the total content of the Universe. If the dark matter 
would be a dark fluid of baryonic matter, such percentage would be quite above the value settled by those constraints. 
Thus, even if the baryonic dark fluid can not be discarded by dynamical methods, the cosmological constraints make 
it unlikely, a fact which might strength the case for exotic type of dark fluids or for other type of dark matter. 



C. Cosmological Constant 

For a cosmological constant, A, = A g pv . We thus have T tt = -e^K, T pp = T zz = e r 2 ^-^ A, T vv = fi 2 A, 
thus from Eq.(f47|), we obtain 

li+(^)) 2 J A=0 ' (61) 

In this way we see that, within our approximation, a non zero cosmological constant can not explain the observed 
behavior because implies that the observed tangential velocity had to be equal to 1, i. e., they should be moving at 
the speed of light. Something similar occurs with the scalar field. 
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D. Scalar Field 



For scalar field <j) with potential, T M „ = <^ jM </>,„ — h g^v 4> a 4>a + g^v V(<f>). We have that, due to the symmetry of our 



space time <fi = <P{Pi z ) the T„„ components are: 

T tt = \e 2{ - 2 ^ (0 p 2 + <\> 2 ) - e 2 ^ V(tf>), (62) 

T PP = \ (0 P 2 - <t> z 2 ) + e~ 2 ^ V(<f>), (63) 

T„ = -\ (0 P 2 - 4> z 2 ) + e- 2 ^ V(<f>), (64) 

T pz = (/> p 4>z, (65) 

T w = -\ e- 2 ^ p 2 (cf> p 2 + 4> 2 ) + e~ 2 ^ p 2 V(4>). (66) 



Inserting these components in Eq.([47|), we obtain that, as in the cosmological constant case 

(\t^) v{ * ) = - (67) 

Again, either the particles move at the speed of light, or the scalar field potential is zero at the equatorial plane. 
When the scalar field potential is zero, that is, we have a massless scalar field, Eq.(|47|) is satisfied, and we have 
to go back to the Einstein's equations as in the vacuum case. Again from Eq.(^), we take the simplest solution 
fj, = p, and, as in the vacuum case, we obtain ip = I In p. The last metric coefficient 7, can be solved in terms of the 
scalar field. Thus, to have a complete solution, it is only left to solve the Klein Gordon equation for the scalar field: 
D 2 S+ ~D6 Dp = 0, which turns out to be the same equation for the metric coefficient ip. However, in this case we do 
not have boundary conditions well defined: The space time is not asymptotically flat; it is not known the form of the 
space time near and at the origin, we are only analyzing the region where the curves are flat; there are no conditions 
for the scalar field at the equatorial plane. What can be concluded at this stage, is that the scalar field does remain 
being a candidate for the dark matter, and thus to contribute with about 25% of the matter of the Universe. 



IV. CONCLUSIONS 



We have found that in a static, axisymmetric space-time, a sufficient and necessary condition in order to have a 
flat profile for the rotational curves in a plane of that space-time is that its metric tensor must have the form given 
by Eq. (|27| ) . This form of the metric must be the one required for galaxies in the region where the rotational curves 
profile of stars is flat. It is important to stress the fact that in the derivation of this expression, only the geometry 
of the space-time was involved, thus, it is independent of the type of matter which generates such a geometry; that 
is, whatever the matter might be, Eq. (p7|) must be the form of the line element at the galactic plane for a static 
axisymmetric space time which can be expressed in the Papapetrou form, Eq. ([!]) with u> = 0, thus, this result is not 
only the general relativistic analog of the Newtonian result for the gravitational Newtonian potential, <j){r) ~ 1/r, but 
it can be used for any type of matter, including those which do not have a clear Newtonian expression, such as the 
scalar field. 

With this idea in mind, we proceeded further using the Einstein's equations, which essentially describe the inter 
dependence of matter-energy and geometry. We had to accept some loose of the generality of our results, in making 
the assumption that the definite relation which we obtained for the metric coefficients at the galactic plane, Eqs.(p6|), 
are also valid in a close by region off the plane, and thus determine a relation for the second derivatives of the metric 
coefficients, Eq.(^6|). Within this approximation, we were able to obtain a constraint equation among the components 
of a general stress energy tensor, T af3 , Eq.©. We tested this expression in four types of stress energy tensors which 
included the traditional types of matter which have been used as candidates for the dark matter in the galactic halos, 
such as the perfect fluid or the cosmological constant. We obtained that for the vacuum case, a cosmic string type of 
matter does generates the observed motion of test particles. Even though the cosmic strings are unlikely objects to 
be the Universe, it was a clear example for the fact that the Newtonian behavior of the density, p(r) ~ 1/r 2 , is not a 
necessary condition for describing the observed motion. We analyzed also the static perfect fluid, and it is interesting 
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that in difference with the Newtonian description, within the general relativistic formulation, we are able to obtain 
conditions for the equation of state of the dark fluid. Even though we were not able to reproduce the well known 
Newtonian result for dust-like fluid, due to our approximation for the second derivatives of the metric coefficients, we 
did showed that a dark fluid with an exotic type of matter, is a candidate for being the dark matter. Furthermore, 
this results represents the Quintessence type of matter at a galactic level. We performed the study in the spherically 
symmetric static space time in appendix [V], and when applied to the perfect fluid static case we did recover the 
Newtonian result, namely that the "dark fluid" could be a well behaved dust like fluid. We further analyzed the 
cosmological constant case, which within our approximation implied that it can not be, and the massive scalar field 
case, which again, within our approximation, turned out that it has to be massless, and the massless scalar field also 
remains being a candidate for the dark matter. 

As we have mentioned above, our results are useful for describing the region where is observed the flat behavior 
of the test particles rotating around the galactic center, it is clearly needed to proceed further in order to be able 
to describe the motion in the complete region, from the center to the exterior. Some preliminary results indicate 
that a combination of perfect fluid with baryonic matter and some of the matter analyzed here, could be in the right 
direction, [l2| . Also, the approximation we made for the second derivatives of the metric coefficients, has to be further 
analyzed. It would be of great help, in order to be able to apply the present description in objets out from the glactic 
plane, to have a sample of the profile of velocities of such objects. In the simpler spherically symmetric case, there 
is no need to make such approximation, and more definite results can be obtained, though with more restrictions on 
the geometry. Besides, there are reasons to belive that the dark matter halo is spherical p^j, thus it might be a good 
approximation the analysis made within this symmetry. We studied the scalar field in this symmetry in pjj . 

Finally, as the best way to study dark matter is through its effects on the dynamics of the visible objetcs, further 
studies along the lines presetned in this work can be performed using gravitational lensing or jets. 

Overall, we consider that the analysis presented in this work, is on the right track in order to determine which is 
the type of matter which constitutes the 90% of the matter in the galaxies. 
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In this appendix we present an analogous derivation for the conditions which the constancy of the tangential veloc- 
ity of circular orbits impose on the metric coefficients for the spherical static case. Furthermore, we also present the 
Einstein equations in this case for a general stress energy tensor. 

It is interesting that, due to the symmetries, in this case we do not have to restrict the analysis to equatorial orbits, 
and that the above mentioned conditions give a closed form for the metric coefficient gu- 

We begin with the line element 

ds 2 = -B{r)dt 2 + A(r)dr 2 + r 2 (d9 2 + sin 2 6dip 2 ) (68) 

The Lagrangian for a test particle reads 

2C = -B(r)i 2 + A(r)r 2 + r 2 {6 2 + sin 2 6 tp 2 ). (69) 

We infer the conserved quantities, the energy E = B(r)i, the (^-momentum L v — r 2 sin 2 9<p, and the total angular 
momentum, L 2 = Lg 2 + (-^g) 2 , with Lg = r 2 9. The radial motion equation can thus be written as: 

r 2 + V{r) = 0, (70) 

with the potential V(r) given by 

v ^=-W) { W)-^- l) - (71) 
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Notice that, due to the spherical symmetry, we do not need to restrict the study to equatorial orbits, this last radial 
motion is valid for any angle 6. For circular stable orbits, we again have the conditions, f = 0, V r = 0, and V rr > 0, 
which imply the following expressions for the energy and total momentum of the particles in such orbits: 



p2 _ mrY 

2B(r) - rB(r) r ' [ ' 



2B(r) - rB(r) r ' 

and for the second derivative of the potential evaluated at the extrema 

rB(r) rr B(r) r /„ _ 2rB(r) r s 

V (r) rr \ extr = 2^ B [ „, . — ' . (74) 

On the other hand, in a similar manner as it was presented in the text, we obtain that the tangential velocity, 
( Wc ) = WFj ((^f ) 2 + sm2 ^(^) 2 )' f° r particles in stable circular orbits is given by: 

<™> 

Thus, imposing the observed condition that this tangential velocity velocity is constant for several radii, this last 
equation can be integrated for the metric coefficient g ft : 

B(r) =B Q r 2 ^\ (76) 

with Bq an integration constant. 

In this way, we again arrive to a theorem, stating that: For a static spherically symmetric spacetime, v c , the tan- 
gential velocity of particles moving in circular stable orbits is radii independent if and only if the gu metric coefficient 
has the form g tt = Bq r 2 . 

Notice that in this case, one of the metric coefficients was completely integrated and the other one, A(r) remains 
arbitrary. Also, as mentioned, the analysis made no suppositions on the plane of motion, so the result is valid for any 
circular stable trajectory. 

Finally, we can construct the Einstein tensor and arrive to the following Einstein equations which give us information 
about the type of matter curving the spacetime in such a way that the motion corresponds to the observed one: 

-(rA(rY + A(r)(A(r)-l)) = 8nT tt , 



Mr) 



2 



2(v c f + l- A(r) 
= 87rJ rr , 



r((v c y + l)A(ry-2v c *A(r) 
2Mr) 2 -SttT,,, (77) 

where ' stands for derivative with respect to r. This study was applied for the scalar field in 

In this appendix we present the generalization of the derivation of the constraint equation among the metric coef- 
ficients, Eq. (|25|) , for the stationary case, where uj =fc 0, described by the line element |l|. 

^From Eq.(^|), we express t, and (p. in terms of E,L, and the metric coefficients as 

* - [{pL 2 e-^ - uj 2 ) E - u L] , (78) 
M 

P 2V 

(p= — (uiE + L). (79) 
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Using these equation in the constraints ones, Eqs.(H) we arrive at: 



2 e~ 2qp (1 - e" 2 ^ E 2 ) + (to E + L) 2 = 0, 



-{e- 2 ^) p E 2 + ( r — ) {loE + L) 2 + ^{ujE + L)lo p E = Q, (80) 



M / p A* 



Solving for E and L, we obtain: 

E = e"> M 



L = ^-^(V^8 VA), (81) 

M 



where 

A = 2e~ A ^ (fi p ~ fiipp)(fi p - 2p,ip p ) - (uj p ) 2 + 



±lo p ^(u p f -A^ p e-^ (fip-v^p), (82) 
B = 2er^( f ip-2^iPp) 2 ~2(Lo p ) 2 . (83) 



For the second derivative of the potential V(p) evaluated at the extreme, we obtain: 

Vpp\ extr = -¥——[{* ^ " 4^ Vv> + — {ujp) 2 ) A 

B [i \ A* / A* M 

-(W P - ^ + 2 (V, P ) 2 +3 f - 4^ Vp)B + 
A* V A 1 / A* 



+ g V^-g) ]. (84 ) 
V A* / A* 

On the other hand, using Eqs.(|79|), and (pi]), in the expression for the angular velocity, Eq. (|l5|) we obtain that: 



_ e 2 ^ VA—B 

n -~ (85) 
P 

where ^4 and Z? are given by Eqs.([S3|). 

As in the static case, Following Chandrasekhar |fl4]| , we rewrite the line element given in Eq.([j]) as: 

li 2 e~ 2 ^ ' J 



^2g-4V_ w 2 ^ I I ^ ^2 e -4,ft ,.2 



e - 2 ^-T)(dp 2 + ^ 2 ), (86) 
thus, in terms of the proper time, dr 2 = —ds 2 , we have that 

dr 2_ M 2 e- 2 ^ ^ 2[1 e 4 V(/x 2 e-^-^) 2 fdp 



M 2 e -4^_ w 2 L ^2 ^ M 2 e -4^_ w 2 

^^- i *-^im\m\ ,87, 



/i 2 I \ dt J \dt 

from which we can write that 



2 -2*0 

/- 1 e o 2 n 2i 

It [1 — u J, 



M 2 e -4V 
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where u° — 4^ is the usual time component of the four velocity, and a definition of the spatial velocity, v 2 , again 
comes out naturally in this way. 



2 _ e 4 ^(/i 2 e- 4 ^-^ 2 ) 2 {dtp 



\dt [i 2 e-^-u? 

2 / , x 2\ 



— - {[?) + UJ )• m 

which is the 3-velocity of a particle measured with respect to an orthonormal reference system, it has components: 

v 2 = v^ 2 + v {p}2 + v {z)2 . (90) 
For the <p— component of the spatial velocity we obtain: 

V M = —[(lje*e- i *-w*)Cl-u], (91) 



and substituting Q from Eq.(85), we finally obtain an expression for the tangential velocity of a test particle in stable 
circular motion: 

fie- 2 ^ VA-luVA^B" 

where A and B are given by Eqs.(|S3|). 

Imposing the condition of constancy for all radii, that is Vp — 0, thus — Vc V \ with a constant, representing 
the value of the velocity, from Eq. (^2|), we finally have that: 

B= {l-v^ )2 )F 2 A, (93) 

where F — ^/i 2 e~ 4 ^ — Jj / (jie~ 2 ^' + iSj . This last expression, represents a constraint among three of the 

metric coefficients, and we can express one of them, say u), in terms of the other two: ip, and fi. In this way, we again 
arrive to an iff condition, namely: The tangential velocity of a test particle moving in a circular equatorial motion 
in an axisymmetric stationary background, has a radii independent magnitude iff the metric coefficients satisfy the 
constraint equation (SF 
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